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A hydrodynamic formulation of the evolution of large-scale structure in the Universe is presented. 
It relies on the spatially coarse-grained description of the dynamical evolution of a many-body 
gravitating system. Because of the assumed irrelevance of short-range ( "collisional" ) interactions, 
the way to tackle the hydrodynamic equations is essentially different from the usual case. The 
main assumption is that the influence of the small scales over the large-scale evolution is weak: this 
idea is implemented in the form of a large-scale expansion for the coarse-grained equations. This 
expansion builds a framework in which to derive in a controlled manner the popular "dust" model 
(as the lowest-order term) and the "adhesion" model (as the first-order correction). It provides a 
clear physical interpretation of the assumptions involved in these models and also the possibility to 
improve over them. 
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I. INTRODUCTION 



The standard model to understand the large-scale features of the matter distribution in the Universe after decoupling 
from radiation can be hardly simpler: a collection of many identical point particles interacting with each other via 
the Newtonian gravitational force in an expanding spatial background |h,E| . Structure arises as a consequence of the 
Q\ • gravitational instability of initially tiny density perturbations. This model neglects relativistic effects, which become 
important only at scales of the order of the horizon and beyond, or when dealing with relativistic velocities. The 
model also excludes nongravitational interactions, which are assumed to be relevant only at small enough scales. 

The general solution to the dynamical evolution of this model is unknown due to the mathematical difficulties. 
N-body simulations, which numerically solve the dynamical equations (see Eqs. (|lj) in the next section), have been 
Q^' very helpful in understanding this evolution. In this work I look for an analytical derivation of some of the relevant 
features in the formation of large-scale structures. One is not usually interested in following the detailed path of each 
Vh ' particle, but rather in some general properties that typically depend on the behavior of a large number of particles 
and which in the end are the kind of data provided by observations, e.g., the smoothed density and velocity fields. 
This naturally leads to the use of a coarse-grained or smoothed description: the evolution of a few, "macroscopic" 
variables is isolated by making suitable approximations for the evolution of the neglected degrees of freedom and for 
their influence on the relevant ones. 

The idea of a spatial coarse-graining has a long history in the field of statistical mechanics, where it has proven 
quite successful as a powerful tool for extracting information about the dynamical evolution of many-body systems. 
In this work I explore its application in the context of cosmological structure formation. Although the procedure of 
smoothing a field is widely used in cosmology, this application has mainly had a descriptive purpose. Unlike this, 
I study systematically the dynamical evolution of the smoothed fields. As we shall see, this method has the merit 
of providing a common framework for different models (dust and adhesion) of structure formation: it offers a clear 
explanation of the approximations involved in each model and of their physical meaning, so that it also opens the 
way to systematically relax them and obtain improved models. Starting from the microscopic equations of motion of 
the particles, I derive an exact set of hydrodynamic-like equations for the (coarse-grained) mass-density and velocity 
fields (Sec. O). These equations do not form an autonomous set, but rather constitute the first ones of an infinite 



hierarchy that must be truncated to become a useful tool. I discuss in Sec. Ill how this can be achieved by resort to 



a large-scale expansion, based on the assumption that the dynamical influence of the small scales over the large-scale 



evolution is weak. The lowest-order term of the expansion yields the dust model (Sec. IV), and it corresponds to a 



complete neglect of the structure below the coarsening scale. This leads to an eventual failure of the model (in the 
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form of pancake- like singularities) , which can be prevented by considering the first-order correction in the expansion 
(Sec. [VJ) : it accounts for the dynamical influence of the structure below the coarsening scale, and I show by means of 
boundary-layer techniques that it gives rise to a robust "adhesive" behavior of the same kind as the adhesion model. 
I end up in Sec. [Vj with a discussion of the results. 

II. COARSE-GRAINING THE BASIC EQUATIONS 

The basic model is a system of nonrelativistic, identical point particles which (i) are assumed to interact with 
each other via gravity only; (ii) look homogeneously distributed on sufficiently large scales, so that the evolution 
corresponds to an expanding Friedmann-Lemaitre cosmological background; and (iii) deviations to homogeneity are 
relevant only on scales small enough that a Newtonian approximation is valid to follow their evolution. Let a(t) 
denote the expansion factor of the Friedmann-Lemaitre cosmological background, H(t) — a/a the associated Hubble 
function, and Qb(t) the homogeneous (background) density on large scales, Xj is the comoving spatial coordinate of 
the i-th particle, u^ is its peculiar velocity, and m its mass. In terms of these variables the evolution is described by 
the following set of equations pF (Vj denotes a partial derivative with respect to Xj): 

1 . x 

x, = -Uj, (la) 

a 

ii i =w i - Hui, (lb) 



V,- • w, = — AnGa 



— X^( x * ~ X J') ~ gb 

3T* 



(lc) 



V,xwj = 0, (Id) 

where Wj is the peculiar gravitational acceleration acting on the z-th particle. Finally, Eqs. (nh must be subjected 
to periodic boundary conditions in order to yield a Newtonian description consistent with the Friedmann-Lemaitre 
solution at large scales ||. 

To implement the idea of a spatial coarse-graining one employs a smoothing window W(z): this smoothing window 
should define a bounded region of space and whose inner structure is smoothed out. In App. W I discuss the general 
properties I will require from a smoothing window. Let L denote the (comoving) coarse-graining length scale. A 
microscopic mass density field and a coarse-grained mass density field can be defined respectively as follows: 

g mlc (x,<) = -^3 J2 <5 (3) (x-x,(i)), (2a) 

g(x, t;L) = J^W O^f^j Qmic(y, t). (2b) 

The physical interpretation of the field p(x; L) follows straightforwardly from the properties of the smoothing window: 
it is proportional to the number of particles contained within the coarsening cell of size ~ L centered at x. A 
microscopic peculiar-momentum density field and the corresponding coarse-grained field can be defined in the same 
way: 

j m *c(x,t) = ^3^u i (t)^ 3 )(x-x i (t)), (3a) 

j(x,t;I) = ||Ff^jj mic (y,{). (3b) 

One can introduce peculiar velocity fields \i m ic and u by definition as j = gu and similarly for u m i C . The physical 
meaning of u(x; L) is also simple: it is the center-of-mass peculiar velocity of the subsystem defined by the particles 



inside the coarsening cell of size « L centered at x. Notice that u is not obtained by coarse-graining u m i C : from a 
dynamical point of view, it is more natural to coarse-grain the momentum rather than the velocity, since the former 
is an additive quantity for a system of particles. Finally, one can define peculiar gravitational acceleration fields w m i C 
and w through a coarse-graining of the force: 

QmicWmic (X, t) = — — - ^ W i (*) S ^ ( X - X * (*) ) > ( 4a ) 

^ ' i 

pw(x, t;L) = JjjsW O^J^-j ft»-w mlc (y, t). (4b) 

The field w(x) has the physical meaning of the center-of-mass peculiar gravitational acceleration of the subsystem 
defined by the coarsening cell at x. 

From these definitions and Eqs. pk) and Mb), it is straightforward to derive the evolution equations obeyed by the 
coarse-grained fields g and u (from now on, d/dt is taken at constant x, and V means partial derivative with respect 
to x): 

!!+3i^=-iv-(eu), (5a) 

^+4Hgu = gw--V-(guu + n), (5b) 

at a 

where a new second-rank tensor field has been defined: 

II(x, t;L)= /|f (^7; ) Qmicij ' *) [u?Mc(y ' *) ~ u(x ' t] L )][ u mic(y, t) - u(x, t; L)]. (6) 

The field II(x) is due to the velocity dispersion, i.e., to the fact that the particles in the coarsening cell have in general 
a velocity different from that of the center of mass. The trace of II is proportional to the internal kinetic energy of 
the coarsening cell, that is, the total kinetic energy of the particles in the reference frame of the center of mass. 

The physical meaning of Eqs. (Ha) and (pp) is simple: they are just balance equations, stating mass conservation and 
momentum conservation, respectively. The term V • II represents a kinetic pressure due to the exchange of particles 
between neighboring coarsening cells (just like in the ideal gas) and it has the same physical origin as the convective 
term V • (guu), i.e., a nonlinear mode- mode coupling of the velocity field. The difference is that the convective term 
couples only modes on scales > L, while the velocity dispersion term codifies the dynamical effect of the coupling of 
the modes on scales > L with the modes on scales < L. The term gw codifies the gravitational interaction between 
the coarsening cells and it is shown later that it can be split in a similar manner into a contribution due to the large 
scales and another due to the coupling of the large scales with the small ones. Although Eqs. <M) look similar to the 
ordinary hydrodynamic equations, there is the important difference that these equations are exact: as one changes the 
smoothing length, the fields g, u, w, II change but in such a way that the equations remain the same (for example, 
upon increase of the smoothing length, part of the dynamical effect described by V-(guu) is shifted towards V • II). 
This property is reflected in that the equations are not an autonomous system for g and u. In fact, they are just the 
first ones of an infinite hierarchy, as can be checked by computing the evolution equations for the fields w and II (see 



Eq. (Bl), for example). To obtain a useful set of equations, it is necesary to truncate this hierarchy by looking for a 



functional dependence of w and II on g and u. This will be the task of the next section. 

III. THE LARGE-SCALE EXPANSION 

In this section I discuss the closure of the hydrodynamic hierarchy at the level of Eqs. (J5|). When the particle 
interaction is dominated by a fast-decaying ( "collisional" ) force, as in normal fluids, this truncation is achieved by the 
assumption of local equilibrium (see, e.g., Refs. |4]||]). In this case, the interaction determines a privileged smoothing 
scale L and it drives the coarsening cells of size L towards an approximate internal thermal equilibrium, as if isolated 
from each other. The evolution of the large scales (3> L) is then ruled by the interaction between neighboring 
coarsening cells. In this way, for example, Eq. (gb) becomes Navier-Stokes' equation. 

One cannot, however, apply this approach to the problem in hand: a system with a long-ranged, unshielded 
interaction such as gravity does not obey the usual thermodynamics (see, e.g., the brief review in Ref. || and the 



more technical remarks in Ref. jj|). Moreover, this long range implies that the interaction is self-similar and does 
not pick up itself a favored coarsening length. One must therefore make use of a different approach to close the 
hydrodynamic hierarchy. 

For this purpose, I introduce the large-scale expansion, which relies on the assumption that, in the context of 
cosmological structure formation, the evolution of the large scales is weakly influenced by what is going on in the 



small scales. This assumption is further discussed in Sec. VI. Here I show how to formulate this idea in order to write 



the helds w and II in Eqs. (|5|) in terms of g and u. Let a tilde denote the Fourier transform of any field: 

0(k)= / dxe lk - x 0(x), (7a) 



v 



^(x) = I^e- k ' x 0(k), (7b) 

k 

where V denotes the volume of periodicity for Eqs. (|l|). Then the definition pp) can be written as g(k; L) = 
W(L\i.) £> m j c (k). Making use of the Taylor-expansion ( [A3] ) and formally inverting it, one gets: 

£ mic (k) = [1 + ±B (L kf + o(L kf] £(k; L), (8a) 

g mic {*) = [^-\b (I V) 2 + o(L V) 4 ] e (x; L). (8b) 

Consider first the field w. Fourier-transforming its definition (0) and using Eqs. (pic) and (pH), one finds: 

gw (k; L) = — ^2 — Pmic(q) £w(k - q) W(L k). (9) 

q#0 q 

This expression collects the mode- mode coupling via gravity of the scales ~ k~ l with every other scale. Take now a 
large scale k~ x 3> L: the assumption of weak coupling to the small scales then means that in the summation in Eq. 
(0), the main contrib utio n arises also from the large scales, q^ 1 , |k — q| _1 S> L. One is then justified to insert the 
expansions (pk) and (A3) into Eq. (M). Transforming back into real space yields the result: 



gw = pw m/ + BL 2 (Ve'V)w m/ +o(LV) 4 , (10a) 

V-w m/ = -4irGa(g-g b ), (10b) 

V x w m/ = 0. (10c) 

Thus, the unknown field w has been written as a functional of the coarse-grained density field. The combination L V 
(<-> L k) can be formally viewed as a parameter which measures the influence of the small scales over the large scales 
and which has been assumed to be small: hence the name large-scale expansion. The field w m ^, which can be called 
the macroscopic gravitational field, represents the gravitational field created by the monopole moment of the matter 
distribution in the coarsening cells, i.e., as if they had no spatial extension (L — 0). It obeys the equations that one 
would have naively guessed and we see that this is not the whole story: there exist a correction due to the higher 
multipolc moments that can then be properly called tidal correction and it represents the coupling of the large to the 
small scales induced by gravity. Hence, this decomposition into a macroscopic field and a tidal correction is analogous 
to the decomposition in Eq. (||b) into convection, guu, and velocity dispersion, n. 

The same procedure can be now applied to the field n. Now, j m j C = Qmic&mic obeys an expansion like (jqb), but 
VL m ic does not, because it is not defined by a straightforward coarse-graining. What one has in this case is slightly 
more involved: 



_ irmc _ [I - \B (LVf + 0{LVf]] 
U " MC grmc [l-iB(LV) 2 + (LV) 4 ]p 



-B (IV) 2 - BL 2 ^ -V + o(LV) 4 

2 g 



(11) 



Because of this minor complication, it is simpler in this case to work directly with the representation in real space, Eq. 
(0) , into which the expansion (nw and the one corresponding to j m i C are now inserted. Because the smoothing window 



effectively restricts |x — y| < L, the fields within the integral evaluated at point y can be consistently Taylor-expanded 
around point x to order [L V) 4 . The final result reads: 

n = BL 2 g(3 ! u)(a i u)+ (LV) 4 , (12) 

where a summation over the index i is implied. Taking this result and ( |10| ) into the hydrodynamic hierarchy (H), and 
dropping the assumed small corrections o(L V) 4 , I finally obtain 

^+3Hg=--V-(gu), (13a) 

at a 

^^- + 4Hgu = gw mf - -V ■ (guu) + BL 2 \(Vg- V)w m/ - -V • [^u)(^u)] 1 , (13b) 

V-w m/ = -4%Ga(g-g b ), (13c) 

V x w m/ = 0. (13d) 

This is an autonomous system of equations for the two coarse-grained fields g and u. Compared to the hydrodynamic 
equations of a normal fluid, we see that, because of the long range of the interaction, its contribution pw to the 
equation for momentum conservation cannot be written as the divergence of a tensor, i.e., as a nonideal correction to 
the kinetic contribution represented by IT. Also the expression to lowest-order for II is another evident difference. In 
the next sections, I explore the dynamical evolution described by Eqs. (J13J) . 

IV. THE DUST MODEL 

In this section I consider the lowest order in the large-scale expansion. This corresponds to formally setting L = 
in Eq. (|l3"|b), thus yielding: 

^+3Hg+-V-(gu)=Q, (14a) 

at a 



du 

~dt 



Hu+-(u-V)u = w mf , (14b) 
a 

w m/ = -4iTGa(g-g b ), (14c) 

V x w m/ = 0. (14d) 



This is the popular and throughly studied dust model (see, e.g., Refs. |l|,|2|]7| ) . The large-scale expansion provides 
a clear picture of the approximations leading from Eqs. ([lj) to this model: by setting L = one assumes that the 
coarsening cells can be thought of as "big particles" lacking completely an internal structure (Fig. fil). This implies 
neglecting (i) the velocity dispersion II compared to the convection term guu; and (ii) the spatial extension of the 
cells and thus the tidal correction compared to the macroscopic gravitational field w™-^ . 

It will be useful to briefly review some results for the dust model. In the linear regime of small fluctuations about 
the homogeneous cosmological background, Eqs. ( |l4| ) can be linearized and the resulting set of linear differential 
equations solved. In particular, in the long-time limit (but still within the linear approximation), the peculiar velocity 
and the gravitational acceleration are related by the condition of parallelism: 

w™/(x,t)=F(i)u(x,i), (15) 

with the function 

F(t) = 47TGg b (t) b J§- > 0, (16) 

b{t) 



where b(t) is the growing mode of the (small) density contrast, i.e., the growing solution of the equation b + 2Hb — 
4TrGg b b = 0. 

The linear evolution eventually breaks down due to the growth of inhomogeneities by gravitational instability. The 
solution to the fully nonlinear equations (|14| ) is not known, but a successful approximation in this regime is the 
Zel'dovich approximation (ZA hereafter) pH3fi, which turns out to be an exact solution of Eqs. ( J14| ) for some highly 
symmetric configurations. In the most general case, it can be understood as the extrapolation of the parallelism 
condition ([l5|) into E n. dl4[ b) J10[ . More properly, considering the way Eqs. ([li]) were derived, the actual approximation 
is the truncated ZA ll|7since the density and velocity fields have been smoothed on a scale L. Comparison with 
N-body simulations [IE shows that the truncated ZA performs substantially better than the original ZA: this is no 
wonder within the present coarse-graining formalism, since the particles that must obey the "dust" evolution (and 
thus be moved according to the ZA) are not the N-body particles of the simulations, which follow Eqs. (pi), but rather 
the "big particles" which represent the coarsening cells. 

If one introduces a rescaled velocity field v = u/ab, and uses b(t) as the new temporal variable, then Eq. (pT|b) 
becomes in the ZA: 

g + (v-V)v = 0, (17a) 

V x v = 0, (17b) 



together with the irrotationality constraint following from Eqs. ( |15| ) and fll4|d). Hence, the problem reduces to the 
curl-free evolution of a fluid under no forcing at all. The solution to this equation is then inserted into the continuity 
equation (|14|a) to yield the evolved density field. As is well-known, however, the ZA has the problem of giving rise 
to singularities in the fields, mainly sheet-like singularities or "pancakes" B at which V • v — > — oo and g — » +oo. 
Beyond this moment, the ZA ceases to be valid. The reason lies in the nonlinear, convective term in Eq. &), which 
deforms the initial velocity field and generically leads to a multivalued velocity field (shell crossing in the cosmological 
literature) . This feature is likely not exclusive to the ZA but a generic property of the fully nonlinear dust model, 
Eqs. (|14|), as checked, e.g., by the application of Lagrangian perturbation techniques (see |Q and refs. therein). The 
generation of singularities by the convective term is not prevented by the gravitational acceleration w m ^, which, on the 
contrary, favors this process. The emergence of singularities signals the unsuitability of the dust model beyond that 
time: in the next section I investigate how the correction to dust following from the large-scale expansion regularizes 
the singularities. 

V. THE ADHESION MODEL 

The singularities of the dust model arise because the approximation of the coarsening cells as "big particles" is 
bad when the density of these "particles" is large (formally infinite): the interaction between them is no longer the 
simple macroscopic gravitational force jw" 1 ', but it becomes more and more dependent on the internal structure of 
the coarsening cells. Therefore, a first step to improve the dust model is to take into account the first-order correction 
following from the large-scale expansion and study Eqs. (|l3|). 

Obviously, the nonlinearities represented by this correction make it even more difficult to solve these equations 
than in the case of the dust model. Hence, in order to understand the effects of the correction on the dust evolution, 
I introduce some further simplifying assumptions. The idea is to consider the limit BL 2 — > + in Eqs. (13h, so 
that the correction is irrelevant everywhere except at those places where the dust evolution would predict a singularity. 
Therefore, the evolution will follow the ZA almost everywhere and the parallelism approximation ( |l5| ) will be good: this 
implies in particular, via Eqs. (jlgc) and (jlSJd), that u is also curl-free and that V • u ex — (g — Qb). Since the correction 
to dust in Eq. (p"3|b) will be effective only near the singularities and these correspond in the ZA predominantly 
to pancakes, the correction can be computed under the approximation of a local plane-parallel collapse: g and u 
change only along the direction n of the local plane-parallel collapse and diUj « (V • ujriirij. Making use of these 
approximations, the correction term can be simplified as follows: 

-{Vg- V)w m/ w -(V-u)Vp«FV 2 u (oc -Vg), (18a) 

(J {J 

— V-[e(ftu)(ftu)]w— V[e(V-u) 2 ]«-(V-u)V 2 u (oc-g\/g), (18b) 

ag ag a 



where it has been taken into account that the correction will be effective only at those places where g ~ (— V-u) — ► +00. 
In this limit, the contribution (jlqb) due to the velocity dispersion will be dominant over that of the tidal correction 
(jlga). Inserting this result in Eq. (|13p) and in terms of the variables v and b previously introduced, I finally arrive 
at the following model: 



~db 



+ (v V)v =z/|V-v| V 2 v (v = WL 2 



H 



(19a) 



V x v = 0, 



(19b) 



This is a closed set of equations for the velocity field: once solved, this field is used to compute the evolved density field 
from the continuity equation ( |13|a ), much in the same way as with the ZA. The coefficient v was called the gravitational 
multistream viscosity in Ref. |12| . Because it multiplies the highest-order derivative in Eq. (J19|a) , one cannot simply 
drop this term in the limit v — > + , but one must rather apply the techniques of the boundary-layer theory to study 
this limit. The mathematical handling is collected in the next subsection; here I simply quote the conclusion that, as 
expected, the correction regularizes the pancake-like singularities predicted by the dust model, which become robust 
structures, where more and more mass gets "adhered" . Another conclusion is that this "adhesive" behavior is rather 
insensitive to the detailed functional dependence on V • v of the factor in front of V 2 v: hence, the same behavior 
arises in particular if one simply drops the |V • v| to obtain a linear correction. But then, one recovers the original 
adhesion model ]MBUl3], which was introduced as a phenomelogical (and, when compared against N-body simulations 
fL3[ , very successful) correction to the ZA to go beyond the epoch of formation of singularities. 

In the framework of the large-scale expansion, the physical explanation of the adhesive behavior is clear (Fig. 
g). The corrections to dust (|l^) behave as a force in the opposite direction to the density gradient, so that in the 
neighborhood of a pancake there is a competition between the inflow of matter driven by the gravitational attraction 
and this "repulsion" , which prevents the formation of a singularity. The dominant contribution as g — > +00 is the 
correction (|l8p) due to the velocity dispersion and represents the conversion of "coherent" streaming kinetic energy 
into "disordered" internal kinetic energy within the coarsening cells |l2|,|l4],|l5[ (the mechanism has the same origin 
as the pressure and viscous-like forces in a gas). The subdominant contribution (|lqa) from the tidal correction 
also exhibits this "repulsive" behavior, and means that the macroscopic field w m f overestimates the gravitational 
attraction between coarsening cells. 



A. Application of the boundary-layer theory 



In this subsection I apply the theory of boundary layers fllq] to the model (fL9f) . The physical picture is that one 
can study the evolution almost everywhere setting v = 0. However, there arise regions ("shocks") spatially well 
separated from each other and of vanishing volume (as v — ► + ) where v has a discontinuity, |Vv| diverges and the 
effect of the term multiplied by v must be taken into account. The fact that these shocks correspond mainly to 
plane-parallel structures greatly simplifies the technical handling, since it reduces the problem to the solution of an 
ordinary differential equation. 

The purpose is to study Eqs. ( |l9[) in the vicinity of a pancake, so that one takes a point at the pancake and 
introduces a new coordinate system moving rigidly with it and with one of the axis normal to the pancake. In this 
new noninertial reference frame Eq. (|19a) must be appended with a term A which collects the acceleration due to the 
inertial forces, and the velocity v is now understood as relative to this new reference frame. New rescaled coordinates 
are defined as x' = x/e and the idea is to take e — > + in such a way that the fields and their derivatives (and 
consequently also the thickness of the pancake) remain finite in terms of the new coordinates x' even in the limit 
v — ► + . (The prime of the rescaled coordinates will be dropped hereafter to simplify the notation). Let n, Xn 
denote the coordinates normal and parallel to the plane of the pancake, respectively, and v n , V11 the corresponding 
components of the velocity v. Let V|| denote the gradient with respect to xy . In terms of these new coordinates, Eqs. 
( |l9| ) now read: 



dv ale 







xn • V, 



d_ 
an 



•Vi 



1— 1 
v = e ' v 



dv n 

dn 



V 






7-2 



dn 2 



+ Vfl I v, (20a) 



d „ 



X (e n V n +V11) =0, 



(20b) 



where I have generalized the model (|l9|) by allowing the prefactor of V 2 v to be an arbitrary function of V • v 
whose asymptotic behavior as V • v — > — oo is characterized by the power 7 (the case 7 = 2 recovers the original 
adhesion model). It has also been considered the case that v (and hence, the rescaling factor e) may depend on time. 
The motivation for this generalization is twofold: (i) This kind of behavior corresponds to the models discussed in 
Refs. 112,04,0], where the hydrodynamic hierarchy (H) is closed by simply neglecting tidal corrections altogether and 
assuming 11™= ng^Sij, k — > + . (ii) The computation of the trace of II from N-body simulations also yields this kind 
of behavior |18| . 

The factor e must be chosen so as to absorb the explicit dependence on v in Eq. (pp|a) and thus to render the 
coefficient of the highest-order derivative of order unity, namely e — i/ 1 /^~ 1 ^> and requiring 7 > 1, so that e — > + . 
This latter constraint states that the correction to dust must grow fast enough with |Vv| so as to succesfully oppose 
the gravitational attraction and overcome the formation of a singularity. Eqs. (EG) are now simplified by keeping only 
the dominant terms in the limit e — » + : first, one can set V|| — » 0, because in terms of the rescaled coordinates the 
pancake looks like an infinite plane and the spatial variations along it take place over an infinitely large lenght scale: 
only the derivative normal to the pancake is relevant. Also, the acceleration terms dv/db and A are finite in the 
neighborhood of pancakes. Finally, de/db is of order e at any finite time. Therefore, Eqs. (|2fj) reduce to the following 
system of ordinary differential equations: 

fin rli) f)^m 

(21a) 



9v n 


dv n 


7- 


' d 2 v n 


dn 


dn 




dn 2 


dvu 
dn 


0. 





(21b) 

The physical interpretation of these equations is straightforward. The tangential velocity, vn , is smooth at the 
pancakes and so constant in the rescaled coordinates. The normal velocity, which looks discontinuous in the physical 
coordinates, is determined in rescaled coordinates by the balance between the convective transport towards the 
pancake (governed in turn by the gravitational attraction of the pancake) and the outwards "pressure" due to the 
velocity dispersion. 

The boundary conditions to be imposed to these equations are: (i) v n (0) = 0, meaning that the rescaled coordinate 
system is centered at the pancake and moves with it; (ii) dv n /dn < for any n, so that there is an inflow of matter 
towards the pancake (and g > 0, see Eq. (23)); (iii) dv n /dn — > as n — > ±00, so that the derivative in physical 
coordinates is not divergent outside of the pancake. Integrating twice Eq. (pl|a) with these boundary conditions, one 
finds an implicit solution v n (n): 
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where V > is an integration constant. Knowning v n (n), one can get an ap pro ximation to the pancake density profile 
by applying the parallelism condition ( |l5[ ) together with Poisson equation (|13]c) in the limit e —* + : 

b_Qb dv n 
e dn 
whose solution reads 



Q = 



(23) 



g(n) 
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n{nf 
V 2 



(24) 



A detailed study of the solution (G3) leads to the general picture shown in Fig. Bt the fields remain univalued at all 
times and no singularity arises. The parameter A — + + measures the pancake thickness in the physical coordinates. 
This study also shows that for the cases 7 > 2 the solution approaches so fast its asymptotic values at large distances 
that in fact dv n /dn = at two finite values of the normal coordinate, n + and n_ = — n + . In such cases, the solution 
( p2| ) must be replaced beyond these points by: v n {n) — (— )V, n < n_ (> n + ). 

From the general solution (p2|), one can compute the behavior of the solution v n {n = z/e) in the limit s — > + for 
a fixed value zq of the physical coordinate z — en: 

v n = V(ilz<z ), 0(i£z = zo), -V(ilz>z ), (25) 

independently of the value of 7 and the temporal dependence of v. The conclusion is therefore that the adhesive-like 
behavior is a robust property of the whole family of models (EG) parametrized by 7: whenever the dust evolution 
(y = 0) predicts a singularity, this must be replaced by the "adhesive prescription" (p5|), which leads, via the continuity 
equation (^3|a), to a steady accretion of mass. 



VI. DISCUSSION 

I have derived and applied a hydro-dynamic- like formulation for the process of cosmological structure formation: this 
is a nontrivial statement, because one could believe that the collisionless nature of the basic model (Q) renders a "fluid" 
description impossible after shell-crossing, and that then one should employ a different approach, e.g. the BBGKY 
hierarchy. I have shown, however, that a "fluid" description is feasible after the breakdown of the dust model: what 
makes a difference with "down-to-Earth" fluids is how the closure of the hydrodynamic hierarchy (|5|) is achieved. For 
this pupose, I have introduced the large-scale expansion, which builds on the assumption that the large-scale evolution 
is insensitive to the small scales. This assumption lies behind many reasonings in the cosmological literature: for 
example, behind the idea that on large scales the evolution follows a Friedmann-Lemaitre model, regardless of the 
small-scale inhomogeneities, and also behind the confidence on cosmological N-body simulations, where each N-body 
particle is so massive that it must correspond in the real world to a full structure in its own, composed of many 
smaller particles. The plausibility of the assumption can be argued on the basis of the long range of gravity: the 
evolution can be expected to be dominated by the large scales provided there is enough large-scale power initially. 
(One can then also expect that the validity of the large-scale expansion should depend on the initial and boundary 
conditions). The good agreement with N-body simulations of the truncated ZA and the adhesion model p 13 1, in 



which small-scale structure is completely disregarded, can be viewed as a support of this large-scale dominance. 

From the derived hydrodynamic equations, I have shown how the dust model arises as the lowest-order term in the 
large-scale expansion, while the first-order correction gives rise to a model which can be reduced to the adhesion model 
by further approximations. This derivation provides a clear physical interpretation of the two models: the "particles" 
which are assumed to follow the dust model are not such but have an internal structure, and the "interaction" between 
these "particles" due to its internal structure explains the adhesive behavior. By applying boundary-layer techniques, 
I found in particular that this behavior is a robust property of the model, in the sense that it arises quite independently 
of the detailed dependence of the velocity dispersion and the tidal correction on the density and velocity fields. 

It is interesting to compare this work withprevious, related works dealing also with a hydrodynamic- like formulation 
of large-scale structure formation. In Refs. JljJT^], the case is studied in which the hydrodynamic hierarchy is truncated 
by a phenomenological ansatz which writes the corrections to the adhesion model as a stochastic term (a noise). The 
adhesion model, in turn, was justified in Refs. [|lJ,[l4|,[l7J (see Ref. p{| f° r a relativistic generalization) by disregarding 
the tidal correction altogether and by letting the velocity dispersion be given as Ily = ng^Sij. The robustness 
of the adhesive property already mentioned explains why this behavior also arises in the models that follow from 
this truncation. It is particularly instructive to compare with the work in Ref. JL4J], where a truncation relying on 
something else than phenomenology was studied. I have been able to close the hydrodynamic hierarchy with less 
restrictive assumptions and thus shown that the adhesive behavior can still be recovered when assumptions (A2) 
(isotropic II) and (A4) and (A5) (further restrictions on the form of II) in Sec. 4 of Ref. pf are dropped. Assumption 
(Al) seems in practice equivalent to the large-scale expansion, while I also employed assumption (A3) (parallelism) 
in order to derive adhesion-like models. But there is also a major improvement compared to that work: I do not 
assume the mean-field approximation from the outset (and thus the Vlasov equation, which was the starting point in 
Ref. p|). In fact, I could derive an expression for the tidal correction and show that it also behaves "adhesively". 
As explained in more detail in App. ||, this correction is also the origin for the discrepancy between expression (|1J) 
and the relationship Ely = ng 5 ^ 3 Sij derived in Ref. p[ . 

Finally, I have shown that the derivation of adhesion-like models requires in principle further assumptions than 
the simple large-scale expansion. Thus, one could improve on these models by relaxing those assumptions and rather 
studying Eqs. ([13). An important difference between the adhesion-like models and Eqs. (O) is that the corrections 
to dust in the latter generate vorticity, even if it is initially absent: this may be a relevant feature when modelling 
galaxy formation. Another difference is that now there is no need in principle to retain the condition BL 2 — > + : 
pancakes and other singularities are no longer of vanishing volume, but have an inner structure whose evolution could 
be studied with Eqs. (|l3|). For such purpose, the role of the length L must be better understood. 
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APPENDIX A: SMOOTHING WINDOWS 

A smoothing window W{z) pjj should behave as a window that defines a bounded region of space and as a smoothing 
filter that erases the structural details inside this region. These conditions are implemented by requiring W{z) and 
its Fourier transform W(q) to decay faster than any power as z — » oo or q — » oo, respectively. The smoothing window 
acts as an integral kernel: if </>(x) is a given field, then the associated coarse-grained field 4>l(x) over the scale L and 
its Fourier transform (f>i (k) are given by 

fe(x) = i|dyff^)^), (Aia) 

0L(k) = W"(ik)^(k). (Alb) 

Therefore, the coarse-grained field at point x is just the (weigthed) addition of the original field over a region of size 
rj L around that point. Two more conditions are also required: (i) W(z = 0) = 1, so that the contribution from the 
neighborhood of the center of the window is unweighted; (ii) W(q = 0) = 1, so that the coarse-grained field has the 
same large-scale structure as the original field. This latter condition implies the following normalization: 

dxW(|x|) = l, (A2) 

so that liniL^o L 3 W(\x.\/L) = S(x), and lim^^o 0l( x ) = 4>{ x )- 

The condition on the decay of W(q) for large q prevents the window W(z) from having sharp borders, which thus 
becomes a "fuzzy" window. Otherwise, the coarse-grained field </>l( x ) could change in a discontinuous manner as the 
center of the window, x, sweeps the system. In the same way, the fast decay of W{z) for large z (which implies that 
W(q) lacks sharp borders in Fourier space) implies the property of spatial locality for the coarse-grained field. This 
also guarantees the existence of the Taylor-expansion of W(q) at q = in the form: 

W(q) = l-^Bq 2 +o(q 4 ), (A3) 

where the constant 

B = - dzz 2 W(z) = y dzz A W{z) (A4) 

is related to the quadrupole moment of the smoothing window. 

The required constraints on the smoothing window exclude a step function (either in real or in Fourier space). A 
useful function which satisfies the constraints is a Gaussian: W(z) = exp(— ttz 2 ) => W(q) = cxp(— q 2 /Air). 

APPENDIX B: THE EVOLUTION EQUATION FOR n 

In this appendix I study the evolution equation for the velocity dispersion and show that ( |l2| ) is a solution of the 
equation. The purpose is to compare with the result in Ref. [|14| for II, also obtained from the evolution equation for 
this tensor. 

Starting from the definition (J5|) and Eqs. (|l|a) and (|l|b), one can obtain the following equation for the temporal 
evolution of the velocity dispersion: 

—^- + hHIiij = — d k (uiJIij) UikdkUj HjkdkUi d k C ijk + Pij, (Bl) 

ot a a a a 

where a summation over repeated indices is implied, and a second-rank and a third-rank tensor fields have been 
defined: 

V(-x,t;L) = / — W ( — - — j £>,m C (y,£)j[u mic (y,£) - u(x,t;L)][w mic (y,t) - w(x,t;L)] + 

+ [w mic (y,£) - w(x,t; L)][u mic (y,i) -u(x,t;L)]\, (B2a) 
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£(x,£;L) = / — W f — - — j Q m ic(y,t)[vL mic (y,t) - u(x,t;L)][u mic (y,t) - u(x,i;X)][u mic (y,i) -u(x, *;£)]. 

(B2b) 

The term £ accounts for the change of velocity dispersion due to the exchange of particles between the coarsening 
cells, while the term V represents the change due to the gravitational interaction. To get a physical picture, just 
consider the equation for the trace of II (the internal kinetic energy) : then C reduces to the equivalent of the kinetic 
contribution to the heat flux in the usual hydrodynamics, while V becomes the power performed by the gravitational 
interaction. 

Compared to Eq. (4c) of [Of, Eq. (Bl) contains the extra term V, and the reason is that this term drops if tidal 



corrections are neglected. Indeed, if one performs a large-scale expansion of V and C in the same way as explained in 



Sec. Ill, one obtains 

V = BL 2 Q[(diu)(diW mf ) + (d 4 w m -0(3 4 u)] + (L V) 4 , (B3a) 

£ = o(LV) 4 . (B3b) 



When these expressions are inserted into Eq. (Bl), one can check that the velocity dispersion tensor II given by (|1J) 
is indeed a solution. 



In Ref. p4] , Eq. (Bl) (with V = 0, as explained) was solved for the trace of II by imposing the constraint of 



shear-free flow: it is then found that TrU = kq 5 ' 3 , where k is determined by the initial velocity dispersion. This 
must be viewed as a solution valid only for early times; when the tidal corr ect ions grow, the term V becomes relevant: 



II eventually forgets its initial condition, due to the term 57111 in Eq. (Bl), and becomes "slaved" to the density 
and velocity fields, as given by Eq. (112). This is a mechanism similar to the "slaving" represented by the parallelism 
condition (^5() in the linearized dust evolution. 
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FIG. 1. The dust model amounts to neglecting altogether the internal structure (velocity dispersion and spatial extension) 
of the coarsening cell: it is approximated by a "big particle" located at the center of the cell, whose mass is that contained 
within it and whose velocity is the center-of-mass velocity. 







Displacement 

FIG. 2. Sketch of the particle trajectories forming a pancake, based on N-body simulations J22J. The horizontal segment 
represents a coarsening cell. The dashed trajectory corresponds to the dust model and its extrapolation (by the ZA) beyond the 
singularity; it is unable to reproduce the stabilization of the pancake. The solid trajectory represents the adhesion model and 
is more realistic. There is no difference initially, but at the pancake the velocity dispersion becomes very large and produces 
an effective "adhesive" forcing that corrects the dust evolution. 
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FIG. 3. Qualitative aspect of the coarse-grained density and velocity fields near a pancake according to Eqs. (E2J) and (b 
A represents a measure of the pancake thickness. 
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